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ABSTRACT 

J  This  paper  treats  travelling  waves  in  a  heterogeneous,  inviscicl,  nou-diffusive  fluid 
bounded  between  two  horizontal  boundaries.  The  fluid  has  two  incompressible  components 
of  different,  but  constant  density  and  is  acted  on  by  gravity.  The  flow  is  steady  when 
viewed  in  a  moving  reference  frame  and  gives  rise  to  a  quasilinear  elliptic  problem  with 
an  eigenvalue  parameter  related  to  the  wave  speed.  The  small  amplitude  solutions  are 
analyzed  using  a  dynamical  systems  approach.  A  center  manifold  reduction  in  combination 
with  a  conserved  quantity  for  the  flow  is  used  to  parametrise  all  ’small’  solutions  of  the 
full  elliptic  system  in  terms  of  solutions  of  an  autonomous  first  order  ordinary  differential 
equation  for  a  principal  component  of  the  wave  amplitude.  The  result  is  a  characterization 
of  all  small  waves,  irrotational  in  each  fluid,  near  the  critical  speed  for  the  system.  They 
are:  solitary  waves;  surges,  connecting  distinct  conjugate  flows  at  extreme  ends  of  the 
channel;  conjugate  flows;  and  periodic  waves.  , 


AMS  (MOS)  Subject  Classifications:  34C35,  35J60,  76B15,  7GC10 
Key  Words:  Internal  wave,  solitary  wave,  surge,  dynamical  system 


*  Supported  by  the  National  Science  Foundation  and  the  Sloan  Foundation. 

**  Supported  by  the  National  Science  Foundation  and  the  Air  Force  OSR. 


SMALL  INTERNAL  WAVES  IN  TWO-FLUID  SYSTEMS 


C.  J.  Amick  and  R.  E.  L.  Turner 


1.  INTRODUCTION 

1.1.  The  physical  and  mathematical  problem 

This  paper  treats  travelling  waves  in  a  heterogeneous,  inviscid,  non-diffusive  fluid 
bounded  between  two  horizontal  boundaries.  The  fluid  has  two  incompressible  components 
of  different,  but  constant  density,  and  the  flow  is  steady  when  viewed  in  a  suitable  reference 
frame,  to  which  we  shall  restrict  attention.  Further  details  about  the  physical  problem  may 
be  found  in  [4], [7], [9],  and  [20]. 

The  vertical  height  of  the  channel  is  normalized  to  be  unity,  and  the  flow  domain  is 
5  =  R  x  I  where  I  =  (— h,  1  —  h),  h  (E  (0, 1)  being  the  specified  depth  of  the  lower  fluid 
in  an  undisturbed  state.  For  a  nontrivial  flow  the  two  components  are  separated  by  an 
interface  T  =  {(:r,F(:r)  :  x  €  R}  C  5,  where  F  is  unknown.  The  density  function  p(x,y ) 
is  normalized  to  equal  unity  in  the  region  S~  where  y  <  Y (x)  and  to  equal  a  prescribed 
p  (E  (0, 1)  in  the  region  S+  where  y  >  Y(a :).  A  pseudo-stream-function  ip  is  defined  in  S± 
by  the  formula  p1/2 q  =  (ipy,  —ip i)i  where  q  =  (?i,  52)  denotes  the  fluid  velocity.  Although 
q  is,  in  general,  discontinuous  across  T,  the  assumption  of  non-diffusivity  ensures  that  it 
is  tangent  to  T  as  approached  from  within  S+  and  S~\  in  particular,  T  is  a  streamline. 
We  assume  there  are  no  stagnation  points;  that  is,  |q|  ^  0,  whence  all  streamlines  go 
from  minus  to  plus  infinity.  The  top  and  bottom  walls,  {y  =  1  —  h}  and  {y  =  —  h},  are 
streamlines. 

The  total-head  pressure  H  is  defined  by 

(1.1)  H  =  p +  ^\Vxp\2  +  pgy, 

where  p  is  the  pressure  and  g  >  0  is  the  gravitational  constant.  A  calculation  from  the 
Euler  equations  ensures  that  Vf7  and  are  proportional.  In  fact  Vif  =  A xpVxp  in  S± 
so  that  q  •  Vf7  =  0  and  H  is  constant  on  each  streamline.  If  we  assume  that  q  approaches 
a  constant  vector  (c,0)  as  x  — >  —00  (or  +00),  then  the  Euler  equations  show  that  H  is 
constant  in  each  fluid  at  infinity.  Since  streamlines  go  from  —00  to  +00,  it  follows  that  H 
takes  constant  values  H±  in  S*.  Since  |q|  7^  0,  it  follows  that  A  ip  =  0  and  so  the  flow  is 
irrotational  in  S±.  This  will  occur,  then,  if 

(1.2)  lim  F(i)  =  0, 

I  —*  —  OO 
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(1-3)  lim  q(x,y)  =  (c,  0) 

X— *■  — OO 

where  c  is  a  constant.  In  the  original  ’laboratory’  coordinates  one  would  see  a  wave  moving 
from  right  to  left  with  speed  c,  into  an  undisturbed  fluid  with  an  interface  having  height 
zero  infinitely  fax  to  the  left.  In  the  moving  frame  the  pressure  is  hydrostatic  at  x  =  — co, 
and  may  be  normalized  by  p(— oo,0)  =  0.  Then  ( 1 . 1 )-( 1 .3)  yield  H  =  H+  =  pc2/ 2  in  S+ 
and  H~  =  c2 /2  in  S~ .  The  pressure  is  to  be  continuous  across  f,  whence  (1.1)  implies 


(1.4)  1/2{|V0+|2  -  |V0~|2}  =  (p  -  1){c2/2  -  gy}  on  T 

where  V0i  denote  the  limiting  values  of  V0  on  T  from  within  S ±. 

This  paper  considers  a  slightly  wider  class  of  flows  than  those  satisfying  (1.2)-(1.3)  in 
that  we  shall  also  allow  periodic  flows.  Since  small  solutions  are  of  interest,  the  function 
tp(x,y)  should  be  close,  in  some  sense,  to  the  trivial  solution 

(1.5)  ’J'(y)  =  c  f  y/p00(s)ds, 

Jo 

where 


(1.6) 


y  <o, 

y  >  o, 


for  a  constant  p  e  (0, 1).  In  particular,  we  shall  restrict  attention  to  flows  for  which 


(1.7)  0(x,  —  h)  —  \f'(  —  h)  =  —ch;  ip(x,  1  —  h)  =  ^*(1  —  h)  =  — cv/p(l  —  h),  x  €  R. 
We  also  want  the  total-head  pressure  to  be  constant  in  the  two  domains  whence 

(1.8)  Vxj  =  A  ip  =  0  in  5+  U  S~ . 


The  mathematical  problem  is  the  following:  given  h  and  p  in  (0,1)  find  unknowns  T,  c, 
and  0  satisfying  (1.7)  and  (1.8)  with  0  =  0  and  (1.4)  holding  on  the  unknown  interface  T. 

We  shall  consider  only  flows  for  which  no  reversal  occurs;  that  is,  q\  >  0,  or,  equiv¬ 
alently,  0 y  >  0.  For  such  flows  one  can  solve  for  y  as  a  function  of  the  spatial  variable  x 
and  the  material  coordinate  0.  The  advantage  of  this  semi-Lagrangian  formulation  is  that 
the  unknown  interface  function  Y  is  now  y(x,0)  evaluated  at  0  =  0.  The  disadvantage  is 
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that  Laplace’s  equation  (1.8)  is  replaced  by  a  singular,  quasi-linear  equation  for  y(x,J4>). 
Let  Y(ip)  denote  the  function  inverse  to  that  in  (1.5)  so  that  Y(^(y))  =  y,  y  €  I.  As  in 
[4]  a  final  change  of  variables  is  performed:  let  aq  =  x,  x2  —  Y(rJ>),  and 

(1.9)  w(x!,  x2)  =  y(xi,V(x2))  -  x2 


so  that  w  represents  the  deviation  of  the  streamline  height  y(x,V>)  from  its  value  Y (ip)  in 
a  trivial  flow.  For  ease  of  notation,  we  subsequently  write  x  for  x2  and  y  for  x2. 

To  describe  the  problem  in  the  new  coordinates  we  define 


(1.10) 


/i(Vu;)  =  T~~ —  and  /2(Vu>)  =  — ^ 

1  +  Wy  1  +  b 


|Vu>| 


' y  x  t  wy  2(1  +  Wy)2' 

The  problem  becomes  the  following:  find  an  eigenvalue  A  =  and  a  function  uj(x,  y), 
continuous  in  T  =  R  x  /,  satisfying 


(1.11) 


r\  r\ 

—  (/l(Vti,))+— (/2(Vw))  =  0  in  ruT+, 


(1.12)  /2(Vui  )  —  Xw  —  p(f2(S7w+)  —  Au>+)  =0  on  y  =  0, 


(1.13)  ty(z,  —  h)  =  iv(x,  1  —  h)  =  0,  x  €  R. 

Here  T~  =  R  x  (  —  h,  0),  T+  =  R  x  (0, 1  —  h ),  and  ±  are  used  to  denote  values  or  limits 
taken  within  T±.  We  ultimately  use  a  standard  weak  formulation  for  a  continuous  w  which 
is  C1  in  T±.  The  theory  in  [4]  shows  that  such  solutions  have  real  analytic  extensions  to 
T±  so  regularity  is  not  an  issue.  The  condition  (1.12)  is  merely  the  statement  that  the 
pressure  is  continuous  across  the  fluid  interface.  Clearly,  for  any  A,  the  trivial  solution 
w  =  0  satisfies  (1.11 )-( 1.13);  we  shall  be  interested  in  nontrivial  solutions. 

One  might  append  to  (1.11)-(1.13)  the  condition  suggested  by  (1.2)-(1.3)  that  w  — *  0 
as  x  — *  — oo  and,  in  addition,  (i)  that  w  — ►  0  as  i  — »  +oo,  giving  a  solitary  wave;  or  (ii) 
that  w  — >0asx— >—  oo  and  approaches  a  nonzero  conjugate  flow  els  i  — *  +oo,  giving  what 
we  shall  call  a  surge;  or  (iii)  that  w  be  periodic  in  x.  Rather  than  any  of  these  conditions 
we  shall  demand  that  w  (and  suitable  derivatives)  be  small  in  T  and  that  A  be  near  the 
critical  value 

<U4>  a'  =  t h(l  +  rh;)- 
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It  will  be  shown  that  such  solutions  satisfy  one  of  (i)-(iii).  After  (A,  w)  has  been  found,  it 
is  straightforward  to  reverse  the  changes  of  variables  to  arrive  at  rp  and  Y (x).  Therefore 
attention  will  be  restricted  to  (1.11)-(1.13)  for  the  rest  of  this  paper. 


1.2.  Methods  and  results 


In  [4]  there  was  shown  to  exist  an  unbounded  branch  C  of  nontrivial  solutions  (A ,w) 
of  (1.11)-(1.13)  bifurcating  from  the  line  of  trivial  solutions  at  the  point  (Ad,0).  The  set 
C  was  considered  in  the  space  R  x  (C0,1(T)  fl  W  '  (T))  and  satisfied  the  following:  (i) 
A  €  (0,  A d),  (ii)  w  is  even  in  x,  (iii)  the  sign  of  w  on  T  and  of  wx  on  (— oo  x  I)  is  the  same 
as  that  of  the  important  parameter 


(1.15) 


J__  p 

h2  (1  -h)2‘ 


In  [4]  it  was  assumed  that  either  e  >  0,  giving  rise  to  a  global  branch  of  waves  of  elevation, 
or  e  <  0,  giving  waves  of  depression. 

Since  the  branch  C  is  unbounded  while  A  is  bounded  (cf.  (i)),  C  contains  a  sequence 
of  solutions,  unbounded  in  C0,1(T)  fl  W  ’  (T)  .  Under  the  assumption  that  they  remain 
bounded  in  C0,1(T)  (which  says,  roughly  speaking,  that  the  corresponding  interfaces  tend 
not  to  develop  vertical  tangents)  suitable  horizontal  translates  of  the  sequence  were  shown 
to  approach  a  surge  (heteroclinic  orbit)  w  satisfying  (1 . 1 1  )-(l .  13).  The  amplitude  of  Tv 
is  proportional  to  |e|  and  the  corresponding  (explicit)  value  A  differs  from  by  0(e2). 
However,  the  existence  of  a  surge  was  dependent  on  the  assumption  of  boundedness  in 
C0,1(T).  Evidence  for  its  existence  can  be  found  in  the  analysis  of  model  equations  done 
in  [10], [12], [18]  and  in  the  numerical  work  of  [19]. 

In  this  paper  we  not  only  show  the  existence  of  small  surges,  but  characterize  all  small 
solutions  when  |A  —  Aj|  is  small.  A  ’dynamical  systems’  approach  is  employed,  whereby  x 
is  treated  as  a  new  ’time’  variable  and  u>(x,  ■)  as  the  ’state’  of  a  system  at  time  x.  One  can 
then  bring  to  bear  a  center  manifold  reduction,  parametrising  all  ’small’  solutions  of  the 
full  elliptic  system  in  terms  of  solutions  of  an  ordinary  differential  equation  for  a  function 
Q  :  R  — *  R.  This  point  of  view  was  developed  by  Kirchgassner  starting  in  the  paper  [13] 
and  has  been  considerably  expanded  in  recent  years  (cf.  [3],[14]-[17]  and  their  references). 
Here  it  is  shown  (cf.  (1.19))  that  the  dynamics  of  the  small  solutions  considered  are 
governed  by  an  autonomous  first-order  ordinary  differential  equation  for  Q  which  depends 
on  A.  With  this  equation  in  hand  it  is  a  simple  matter  to  determine  all  solutions  and  their 
dependence  on  parameters. 
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The  plan  of  the  paper  is  as  follows.  In  section  2  the  necessary  mathematical  notation 
is  introduced.  The  problem  (1.11)-(1.13)  is  considered  in  the  form  of  a  linear  operator 
(the  formal  Frechet  derivative  of  (1.1 1)-(  1. 13)  at  (A^O))  acting  on  w,  equal  to  a  nonlinear 
function  of  (\,w).  The  linear  operator  is  studied  in  section  3  and  the  nonlinear  problem 
solved  in  section  4  with  the  aid  of  the  contraction  mapping  principle.  The  main  idea  is 
to  seek  w(x,y)  in  the  form  Q(x)t(y)  4-  f?(x,y)  where  t  is  an  eigenfunction  of  a  suitable 
boundary  value  problem  in  the  variable  y  and  the  ’remainder’  R  is  orthogonal  to  t(y)  at 
each  x  with  respect  to  the  weight  p^.  Ultimately  R  will  be  of  higher  order  than  linear 
in  Q  and  A —  A.  An  important  reduction  theorem  (Theorem  4.6)  states  that  R  and  its 
partial  derivatives  are  pointwise  functions  of  Q  and  Q' .  Thus  R(x,  •)  =  Fi(Q(x ),  Q'(x),  A). 
To  obtain  this  striking  reduction  we  are  led  to  pose  the  problem  in  such  a  way  that  the 
natural  space  consists  of  functions  which  may  grow  exponentially  at  infinity.  While  this 
may  appear  rather  far  from  our  goal  of  obtaining  small,  bounded  solutions,  the  success  of 
the  contraction  mapping  hinges  on  working  in  a  suitable  closed  subset  of  these  growing 
functions.  Section  3  is  necessarily  long  due  to  the  technicalities  of  working  with  possibly 
unbounded  functions. 

In  the  nonlinear  equation  satisfied  by  w  one  then  expresses  R  and  its  derivatives  in 
terms  of  Q  and  Q'  to  arrive  at  an  ordinary  differential  equation  for  Q: 

(1.16)  Q'\x)  =  F(Q(x ),  Q'(x),  A),  x  e  R 

where  F  is  even  in  its  second  variable.  Analogous  equations  arise  for  water  waves  in  the 
presence  or  absence  of  surface  tension  [3].  Although  (1.16)  may  be  analyzed  by  scaling 
and  using  the  implicit  function  theorem  (cf.  [5]),  there  is  a  first  integral  of  (1.16)  which 
simplifies  the  analysis,  carried  out  in  section  5.  In  an  important  paper  [7]  on  internal 
waves,  Benjamin  considered  the  flow  force  defined  by 

(1.17)  ^{#  (*, yS>  +  \^^x' y )  “  ^*( x ’  v))  ”  y)y)dy 

for  each  x  €  R,  and  showed  it  was  independent  of  x.  In  terms  of  w  the  result  becomes 

r  U)^  -  yj  2 

(1.18)  /  pooj - t-dy  +  A {p  -  l)ic2(x,  0)  =  constant 

Jl  (1  +  Wy) 

and  when  all  quantities  are  again  expressed  in  terms  of  Q  and  Q1,  (1.18)  yields  a  relation 

(1.19)  Q'(x)2  =  h(Q(x),  A),  x  e  R. 
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Since  in  (1.16)  the  dependence  of  F  on  Q'  is  even,  (1.19)  shows  that  Q"  is  a  function  of  Q 
and  A  alone,  and  (1.19)  is  just  the  first  integral.  Its  analysis  is  straighforward  and  yields 
the  totality  of  nontrivial  small  waves:  solitary  waves,  surges,  conjugate  flows,  and  periodic 
waves.  For  the  problem  of  solitary  water  waves  considered  in  [3]  equations  analogous  to 
(1.16)  and  (1.19)  hold  when  the  Bond  number  vanishes  or  exceeds  1/3.  In  the  range  (0, 1/3) 
(1-16)  is  replaced  by  a  fourth-order  equation  and  (1.19),  by  a  relation  among  third  and 
lower  order  derivatives. 


2.  NOTATION 


Let  I  denote  a  bounded  open  set  in  Rn  or  a  point  and  let  T  =  R  x  J.  Thus  when  I  is  a 
point  we  merely  identify  T  with  R.  Let  CJ(T)  denote  the  space  of  continuous  real- valued 
functions  on  T  with  continuous  derivatives  through  order  j.  For  h  in  C°(T)  and  0  <  a  <  1 
we  measure  the  Holder  constant  of  h  ’at  x'  by 


(2.1) 


(. Hah)(x )  =  sup 


I h(x,y)  -  h{x\y'\ 


l(*. »)-(*'. »'l<i  \(xiU)  {x  iV  I 

ve/ 


where  |  |  denotes  the  Euclidean  distance  in  T.  For  each  real  number  p  we  define  a  norm 

(2.2)  SUP  e-'‘|l||a',ft(.T,!/)|  +  V  sup 

W<i‘€R',e'  l0l-l*eR 

using  the  standard  multi-index  notation  for  derivatives,  and  let 

(2.3)  =  (fc  €  C>  :  <  oo)}. 

Let  C/’“(T,  Rn)  denote  Revalued  functions,  normed  by  the  supremum  of  the  norms  of 
the  components;  and  C/,  the  space  with  the  HQ  estimate  omitted  from  (2.2).  For  6  >  0 
let 


(2.4)  =  {u  €  C£Q(T)  :  \d*u\  <  6, \f3\  =  k,  1  <  k  <  j;  Had0u  <  8 ,  \0\  =  j} 

and  Cjj(T,  Rn  ),  the  vector- valued  version,  or  simply  when  the  context  is  clear.  Note 
that  CJ,ac  is  a  closed  subset  of  C}'a  and  that  u  £  C]’°c  can  be  unbounded. 

H,o  t* 

For  the  internal  wave  problem  I  will  denote  the  interval  (— h,  1  —  h).  We  shall  repeat¬ 
edly  refer  to  sets  and  functions  restricted  to  the  regions  {y  >  0}  and  {y  <  0}  and  will  use 
the  symbols  +  and  —  for  this  purpoese.  Thus  I+  =  (0, 1  —  h ),  I~  =  (  —  h,  0),  T+  =  R  x  I+ , 


6 


a-fld  T  —  R  x  I  .  For  a  function  u  defined  on  I  or  T,  denote  its  respective  restrictions 
and  if  u±  are  initially  defined,  u  will  denote  the  corresponding  function  defined  on  the 
union  of  the  domains.  In  general,  u  will  not  extend  continuously  to  {y  =  0}. 

Spaces  particular  to  the  internal  wave  problem  are: 

(2.5)  Xl  =  {«<  €  C\T) :  u(x, - h )  =  u(x,  1  -  h)  =  0,  *  €  R), 

the  norm  being  the  sum  of  those  on 

(2.6)  U,  =  x  Cyl\K)  x  Xj, 

and 


(2.7)  Am  =  C°//2(T+)  x  C°/'2(T+)  x  x  C°/'2(T~)  x  C“’^2(R). 

Finally,  denotes  a  product  analogous  to  (2.6)  but  with  subscripts  y.,6  on  each  factor. 

3.  A  LINEAR  PROBLEM 

In  this  section  we  shall  consider  the  linearization  of  ( 1 . 1 1 )-( 1 . 13)  about  (\,w)  = 
(A^,  0).  Formally,  the  problem  is 

(3.1)  Aw±  =  div(gf ,gf)  in  T±, 

(3.2)  p(wy+ -  \dw+)  -  (wy~  -  \dw~)  =  g  on  y  =  0, 

(3.3)  w~(x,  -h)  =  w+(x,  1  -  h)  =  w~(x,0)  -  w+(x,0)  =  0,  x  £  R. 

This  problem  will  be  analyzed  in  suitable  spaces  and  the  results  applied  to  the  nonlinear 
problem  in  the  next  section. 
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3.1.  A  weak  formulation  of  the  linear  problem 
Define 

(3.4)  <  >1=  JpooivW'ytydy  -  \d(l  -  p)x!>(0)cl>(0),  V0,  <jf>  G  C0°°(/) 

and 


(3.5) 

<  -0,^  >r= 


J  J  Pco(y)V^V(f>dxdy  -  Ad(l  -  p)  J  0(x,  O)0(x,  0)dx,  V0, 0  e  C£°(T) 


where  p <*,  is  given  by  (1.6).  Fix  y  >  0  and  G  =  (<7^,5a  ,9% ,92  > 9 )  €  AM.  By  a  weak 
solution  of  (3.1)-(3.3)  is  meant  a  function  €  X*  satisfying 


(3.6) 


<  w,<j)  >t=  [G,  0],  V0eCo°°(T) 


where 


(3.7) 


^=7/ 


Poc(y)(gu92)  ■  V0  + 
T 


{pgt(x,  0) 


-  p2  (x,  0)  -  #(x)}0(x,  0)dx. 


3.2.  The  kernel  of  the  linear  problem 

First  we  consider  G  =  0  in  (3.6).  We  shall  prove  that  if  y  >  0  is  sufficiently  small, 
then  w  belongs  to  a  two-dimensional  subspace  of  X*.  The  dimension  of  this  subspace 
determines  the  order  of  the  differential  equation  (cf.  (1.16))  which  governs  small  solutions. 
A  simple  argument  (cf.  [4])  shows  that 


(3  8) 


max  Aj(l  —  p)r2(0) 
r€W  J^ZFT2 


Jrl  2 

where  W  =  W  ’  (/),  with  the  maximum  attained  by 
(3.9) 


(  x  _  f  My  +  h),  y  <  0, 

'J  ’  \ -mh(l  -  h)~l(y  -f  h  -  1), 


y  >  0, 


and  m  is  chosen  so  that  fj  poo(y)t2(y)dy  =  1.  Further 

Arf(l-p)r2(0) 


(3.10) 


max 


T  €  W,  fjpootT  =  0  fjpoo(r'y 


<  1 
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with  so  restricted.  From  (3.8)  one  sees  that  neither  (3.4)  nor  (3.5)  is  coercive,  while 


(3.11) 


<  t,  $  >=<  r,  <A  >,  +  J 


v  1  2 

is  so  and  extends  to  an  inner  product  on  W  '  (/)  equivalent  to  the  usual  one.  Consider 

2 

the  problem  of  finding  all  nontrivial  solutions  (u,  r)  ERxll’  ’  (I)  of 

(3.12)  <T,4>  >=(!/  +  yj  pooTt,  V<t>  e  W1,2(I). 

Note  that  from  (3.8)  and  (3.11),  v  >  0  and  u\  =  0,  r  =  t  =  <i  is  a  solution.  Since 
the  right-hand  side  of  (3.12)  is  symmetric  and  compact  with  respect  to  <  •,  •  >  standard 
theory  ensures  that  (3.12)  has  a  sequence  of  simple  characteristic  values  0  =  u\  <  v2  <  ^3... 
and  corresponding  eigenfunctions  which  form  an  orthonormal  basis  for  L2 (I)  with 

respect  to  the  measure  poody. 

Assume  w  6  X *  satisfies  (3.6)  with  G  —  0  and  set  <f>(x,y )  =  $(x)tk(y),  where  $  € 
C£°(R).  Although  <j)  ^  Cq°(T),  a  density  argument  shows  it  is  a  suitable  test  function. 
The  use  of  (3.12)  yields 


(3.13) 


/oo  /*oo 

$'(r)A'A.(x)(ii  =  — i/jt  /  $(x)Ajt(x)dx, 
-00  •/  —00 


where  Ak(x)  denotes  the  k-th  Fourier  coefficient  of  w(x,  ■)  : 

(3.14)  Ak(x)  =  J  Poc{y)w(x,y)tk(y)dy,  k  =  1,2,. 


It  follows  that  j4jf(jr)  =  i/kAk(x).  Since  jc  € 

|A*(z)|  <  const. 1  €  R. 


Hence,  if  p  is  restricted  to  lie  in  the  interval  (0,  ■y/i'J),  then  Afc  =  0  for  k  >  2,  and 
Ai(x)  =  £1  -f  f°r  some  pair  =  (£1,(2)  6  R2.  Since  the  form  a  basis,  we  conclude 
the  following: 

Theorem  3.1.  Suppose  p  €  (0,  yF^).  If  w  X *  satisfies  (3.6)  with  G  =  (y^,  gf,  g)  =  0 , 

then 

w(x,y)  =  (<fi  +  £2x)t(y)  inT 
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for  some  £  =  (£1,^2)  €  R2. 


3.3.  The  inhomogeneous  linear  problem 

Next  (3.6)  will  be  solved  for  w  6  X *  in  terms  of  G  (E  AM.  The  parameter  p  will 

eventually  be  restricted  to  lie  in  (0 ,JI)  where  p  =  The  function  w  is  sought 

in  the  form  Q(x)t( y)  +  R(x,y)  where  R(x,-)  is  p ^  orthogonal  to  t(y)  for  each  x  and 

Q( 0)  =  Q'( 0)  =  0.  According  to  Theorem  3.1  the  general  solution  can  then  be  obtained 

by  adding  an  arbitrary  linear  function  £1  +  £2^  to  Q. 

•  •  •  0  1  2 

We  begin  with  solutions  in  the  space  W  ’  (T).  Let  E  denote  the  closed  subspace  of 
0  1  2 

W  ’  (T)  consisting  of  functions  b  satisfying 

Poo(y)t(y)Kxi  y)dv  =  Vx  e  R 

Note  that  E  fl  C£°(T)  is  dense  in  E.  Using  the  preceding  analysis  one  verifies  that  for 
a  €  W1,2(R)  and  b  £  E 

(3.16)  <  a(x)t(y),  <j)  >t=  0,  6  E 

and 

(3.17)  <  b ,  $(x)t(y)  >T=  0,  €  Wj’2(R) 

0  1  2 

so  that  each  w  €  W  1  (T)  has  a  splitting  w  =  a(x)t(y)  4*  b(x,y)  orthogonal  with  respect 
to  <  •,  •  >7’  . 

Let  70  be  in  C^°(R,  [0, 1]),  supported  on  [—2,2],  equal  to  1  on  [—1, 1],  and  such  that 
7o(x)  — 1/2  is  odd  about  x  =  3/2  on  [1,2].  Then  7it(r)  =  7o(x  —  3 k),  for  k  €  Z  is  a  partition 
of  unity  on  R.  Consider  equation  (3.6)  with  G  replaced  by  7 kG  =  (7kdt  ilkd^lkd)-  It 
will  be  shown  to  have  a  solution  of  the  form  ak(x)t(y)  -f  bk(x,  y)  where  bk  €  E.  Then  in 
Theorem  3.2  it  is  shown  that 

OO 

w  =  ^2  +  bk) 

k=  —  o 0 

belongs  to  X *  and  is  the  desired  solution  of  (3.6). 

To  find  the  component  bk  one  solves 

(3.18)  <  bk,  <}>  >t=  [ikG-,  4>],  V<j>€E. 


(3.15) 


J, 
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From  the  inequality  (3.10)  it  follows  that  <•,•>/  is  coercive  when  restricted  to  the  p0 0 
orthogonal  complement  of  t(y)  and,  by  integration,  <  •,  •  >t  is  coercive  on  E.  Since  7* 
has  compact  support,  the  right-hand  side  of  (3.18)  is  easily  seen  to  define  a  bounded  map 
on  E  and  the  Lax-Milgram  theorem  provides  a  unique  solution  bk  £  E  for  (3.18). 

To  determine  ak  one  may  justify  setting  (f>  =  $(x)t(y),  $  £  Cq°( R)  in  (3.6)  where  G 
is  replaced  by  ykG.  The  use  of  (3.17)  and  the  identity  <  t}  t  >/=  0  yield 

/oo  poo  poo  p 

a'&  =  $t(0)jk{pg£(-,0) -£2"(->0) -9}  -  /  /  Poo7k(9u92)  •  ($'<,$*')• 

-00  J —00  J —00  J 1 

The  unique  solution  a*  with  a*(0)  =  0^(0)  =  0  is 
(3.19) 
ajfc(z)  =  - 

+  /  7*(s)  /  Pocgi(s,y)t(y)dyds  -  xyk(0)  I  poogi{0,y)t(y)dy. 

Jo  Jl  JI 


f  (x -■s)7k(s){t(0)(pg2(*,0) -ff3  (s,0)  -  g(s))  +  f  Poo92(a,y)t'(y)dy}ds 
7-00  Jl 


Theorem  3.2.  Suppose  p  £  (0,p)  and  let  wk  =  a.k(x)t(y)  +  bk(x,y)  denote  the  solution 
of  (3.6)  with  G  replaced  by  ykG.  Then 

OO 

w  =  wk 

k=~  00 

is  an  element  of  X*  satisfying  (3.6).  Moreover, 

(3.20)  |u;|Xi  <  const. jGjA((. 


Proof.  Fix  C  €  R,  set  V  =  (£+  1/4,  (  +  3/4)  x  I,  and  U  =  (£,  £  +  1)  x  I.  Let  X  denote 
the  natural  restriction  of  X *  to  V  and  for  u  £  X  take  as  a  norm: 

luk  =  lu+l  cll,-(V+)  +  I  U~\clyl/2(V-)- 

Since  G  =  (gf,  g±,  g)  £  A^,  it  is  immediate  from  (3.19)  that  dk(-)t( ■)  £  X  for  every  k. 
Since  the  function  7 k  has  support  on  (3 k  —  2,3k  +  2],  only  a  finite  number  of  the  ak  are 
non-zero  on  [£,  (  +  lj.  Hence 

m 

£  «*(•)«(•) 

fc=  — m 
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converges  in  X  as  m  — >  oo  to  an  element  a(-)t(-),  where  a  is  given  by  (3.19)  with  -yt 
replaced  by  unity.  From  this  representation  of  a  it  is  immediate  that 

(3.21)  <  const. IG|a„ 

where  the  constant  is  independent  of  (  and  G,  though  does  grow  as  (i  approaches  zero. 
The  rest  of  the  proof  is  devoted  to  estimating  6*.  In  order  to  use  elliptic  estimates 

°  1  j 

we  require  an  equation  for  bk  in  which  all  test  functions  <f>  €  W  ’  ( T )  are  used.  By  using 
a  splitting  <f>  =  $(a:)t(y)  +  ip(x,y)  with  rfi  €  E  we  can  appeal  to  (3.16)-(3.18),  and  the 
equations  satisfied  by  a*  and  bk  to  write 

<  bk,  <f>  >T  =<  bk,  ifr  >T 

=  bfkG,  <f>  -  $<] 

(3<22)  =[-ykG,<t>)-<akt,$t>  T 

=  bkG,  t}-  J  JT  Pooa'kt<t>x 

where  a*  is  given  by  (3.19).  One  verifies  that 

/  It  P°°a'kt<i>z  =  I  Jr  Poo  Ky)bk(x)  J  Poo9i(x,  z)t(z)dz]\l<t>x(x,  y) 

+  J  JTPoot(y)bk{x)t(0){pgt(x,0)  -  gi{x,  0)  -  g{x)) 

+  J  Poo92(x,  z)t(z)dz]<f>(x,y). 

Let  k  denote  the  smallest  integer  with  3k  —  1  <  £  <  3A:  -f  3,  so  that  3 (fc  —  1)  +  2  <  £  < 
(  + 1  <  3 (k  +  2)  —  2.  Then  the  support  of  7 *  is  outside  [£,  C  +  1]  unless  k  =  kork  =  k  +  l. 
If  k  ^  k,  k  +  1  and  (f>  is  restricted  to  C£°(U)  then  from  (3.22), 

(3.23)  <  bk,<f>  >t=  0,  V4>  €  Cg°(U). 

Standard  elliptic  estimates  [2]  give 

(3.24)  Ibklt  <  const.  c“'l|<l|6jt|L2(l/),  k  £  k,  k  +  1 

where  the  constant  is  independent  of  y.,  and  k.  When  k  =  k  or  k  +  1  the  regularity 
theory  gives 


(3.25) 


<COnst.  e  ^  ^  I C7°,  1  / 2 ( C/+ )  bk9i  1 C70. i/a(ty—  )] 

+  IS,!co.‘/2([<iC+i])  j 


<const.  e  ',|cl|6fcUl(i/)  +  const.  |G|a„. 
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The  theory  for  (3.24)-(3.25)  differs  from  the  standard  Schauder  estimates  in  two  ways. 
First,  one  obtains  C1,1/2  estimates  rather  than  the  usual  C2,1//2  estimates,  and  second, 
boundary  conditions  are  posed  on  the  line  {y  =  0}  between  two  regions.  The  bounds  are 
obtained  by  applying  the  theory  in  section  8  of  [2]  for  equations  in  variational  form,  to  the 
functions  pb%(x,y)  +  bk(x,-y)  and  b~£(x,y)  -  bk(x,-y)  for  y  >  0.  The  arguments  are 
technical,  but  straightforward. 

The  next  objective  is  to  estimate  bk  in  L2(U).  Restrict  attention  first  to  the  case 
k  <  k.  It  follows  from  (3.22)  that  bk  satisfies  (3.23)  for  all  <f>  £  Co°((3k  +  2,  oo)  x  I  and  the 
regularity  theory  just  cited  gives  bf  €  C°°((3&  +  2,  oo)  x  I*  as  well  as  the  pointwise  decay 

.  •  ®  i  2 

of  bk  and  V&*  to  zero  as  x  — ►  +oo,  since  bk  is  known  to  be  in  W  ’  ( T ).  In  particular,  bk 
satisfies  the  homogeneous  version  of  (3.1)-(3.3)  on  (3k +  2,  oo)  x  I.  Multiply  equation  (3.1) 
for  bk  by  poobk  and  integrate  over  (x,oo)  x  /  to  obtain 


ld_ 

2  dx 


l 


poob\(x,y)dy 


-  J  J  P<x>\^bk(s,  y)\2dyds 
f°° 

+  \d(l  —  p)  J  b2k(s,0)ds,  x>3k  +  2. 


The  use  of  (3.10)  and  (3.12)  yields 

~  JjPoo\^(s,y)\2dy  +  Ad(l  -  p)b\(sy 0)  <  -u2  J^Poob2k(s,y)dy 
for  all  s  E  R.  It  follows  that 


d_ 

dx 


L 


Poob\(x,y)dy 


<-^n 


Poob2k(s,y)dy 


whence 

~  f  Poobl  <  e2F(3^+2-ri  f°°  j  Poohl 

J  x  Jl  J  3fc+2  J I 

for  all  x  >  3k  +  2,  where  Ji  =  \Jv2/2.  In  particular, 


(3.26) 


J  J  Poob2k  <  const. e2M(3fc  c)  J  J  Poo b\. 


The  use  of  4>  =  bk  in  (3.18)  and  the  Poincare  inequality  yield 

f  f  /*3fc+2  ^  2  f 

/  /  poo b\  <  const.  /  {<?(x,0)2 +  fff(x,0)  +  firf(x,0)  +  j  (g\  +  g])dy}dx 

J  Jt  Jzk-2  Jl 

<  const.  e6^\G\2A 


(3.27) 
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and  this  estimate  holds  for  all  k  £  Z.  Its  use  in  (3.26)  gives 


I bk\Lj(U)  <  const.  e^zk  0+#l3|fcl)|G|A/i ,  k  <k. 

A  similar  argument  yields 

\bk\La(U)  <  const.  e(M(-3fc+C)+M3|fe|)|G|A^  fc>  fc  +  1. 

Since  (3.27)  holds  for  k  =  k,  k  +  1  the  unified  estimate: 

IMl2(1/)  <  const.  e(-pl3*-Cl+M3l*l)|GjA^ 

holds  for  all  k  £  Z.  The  sum  over  k  is  tantamount  to  a  convolution  by  a  kernel  with  a 
decay  parameter  /I,  which  is  known  to  take  a  space  with  weight  0  <  <  /7  into  itself.  The 

result  is 

OO 

(3.28)  ]jT  l&fcl L2(t/)  <  const.  e^|G|A(i. 

*=-oo 

It  is  then  immediate  from  (3.24)  and  (3.25)  that  ]T)|jfc|<m  ifc  converges  in  X  to  an  element 
b  with  |6|^  <  const. |G|A>1 ,  so  that  w(x,y)  =  a(x)t(y )  +  6(x,y)  €  AT  satisfies  (3.6)  for  all 
<£  €  C0°°(V).  Since  C  is  arbitrary,  the  theorem  is  proved,  q.e.d. 

An  operator  M  defined  on  a  function  space  over  T  is  called  translation  invariant  if 

M(ur)  =  (Mu)r 


where  for  all  r  €  R  and  (x,y)  £  T ,  ur(x,y )  =  u(x  +  r,  y).  If  7,1-G  is  replaced  by  7*(jr 
in  (3.6),  the  resulting  solution  will  be  ( 5*)r ,  and  these  sum  to  bT.  Hence  the  linear  map 
which  takes  G  €  AM  to  b  £  X*  is  translation  invariant. 

The  next  theorem  summarizes  the  analysis  of  this  section  and  introduces  some  notation 
useful  for  the  nonlinear  problem  in  the  next  section.  Throughout  the  rest  of  the  paper  it 
is  assumed  that  0  <  p  <  JI  =  -y/i/2/2. 

Theorem  3.3.  There  exist  bounded  linear  maps  L\,  T2  :  AM  — ►  G°’^2( R)  and  L3  :  AM  — ► 
X^  such  that,  given  G  £  AM  and  (£j,£2)  £  R2,  the  weak  form  (3.6)  of  (3.1)-(3.3)  has  a 
unique  solution  w(x,y )  =  Q(x)t(y)  +  R(x,y),  with  ( Q,Q',R )  £  U ^  satisfying 


(3.29) 


Q(x )  =  Ci  +  /  Q',  x  e  R, 

Jo 

Q'(x)  =  £2  +  £iG|5  -f  /  T2G,  x  g  R, 

Jo 

R(x,y)  =  (L3G)(x,y),  ( x,y)£T . 
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Here 

(3.30)  ( hG)(x )  =  J  Poo(y)t(y)gi(x,  y)dy, 

(3.31)  (L2G)(x)  =  t(0){g(x)  +  g2(x,0)  -  pg£(x,0)}  -  j  poo{y)t'(y)g2(x,  y)dy, 

£3  is  translation  invariant,  and 

P<x>(y)t(y){L3G)(x,y)dy  =  0,  x  6  R. 

Remark  3.4.  Theorem  3.3  is  stated  for  p  >  0,  and  it  will  be  important  in  section  4  to 
have  results  for  G  €  Ao-  Since  Ao  C  for  p  >  0,  the  solution  w  =  at  +  6  has  at  and  b  in 
X*  for  all  p  €  (0,  Jt).  As  observed,  the  constant  in  (3.21)  grows  as  p  —>  0.  However,  the 
constants  in  (3.25)-(3.28)  are  independent  of  p  so  there  is  the  additional  estimate 

(3.32)  \L3G\xx  <  const. |G|a0. 


4.  THE  NONLINEAR  PROBLEM 

The  nonlinear  equations  (1.11)  and  (1.12)  can  be  expressed  as 

(4.1)  Aw±  =  div(gi(Vw±),  ^(Vw*))  in  T±, 

(4.2)  p(wy+  -  Adu>+)  -  (wy~  -  XdW~)  =  gCVw*,™,  A)  on  y  =  0, 
where 

(4.3) 
and 

(4-4)  ^(Vw±,  w,  A)  =  pg2(Xw+)  -  02(Vu>“)  4-  (A d  -  A)(l  -  p)w. 

The  intention  is  to  show  that  solutions  of  (3.1)-(3.3)  in  X *  for  which  |A j  —  A|  is  small 
and  for  which  w  and  its  derivatives  are  small  uniformly  in  T,  have  the  decompostion 


WXWy  _  W]  [V^l2 

1+U?y’  92  1  +Wy  2(1  +Wy)2' 
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w(x,y )  =  Q(x)t(y)  +  R(x,y),  where  R  and  its  derivatives  are  pointwise  functions  of  Q(x) 
and  Q'(x).  This  ’center-manifold’  reduction  is  then  used  to  derive  a  differential  equation 
for  Q(x)  to  be  analyzed  in  section  5. 

4.1.  A  transformed  problem 

Set  Q i  =  Q,  Q 2  =  where  /?  >  0  is  a  constant  to  be  chosen  in  due  course  and  let 

£  =  (£i,6)  £  R2  denote  the  value  of  (Qi(0),  Q2(0).  Let  ft  =  (Qi,Q2,  R)  and  p  =  \d  — 
Equation  (3.29)  can  be  expressed  as 


(4.5) 


<2i(x)  =  £i  + /?  /  Q2,  x€R, 

Jo 

Q2(x)  =  6  +  0-%G|5  +  r1  r  L2G, 

Jo 

R(x,  y )  =  ( L3G)(x ,  y),  (x,  y)  G  T 


x  e  R, 


or  in  symbolic  form  ft  =  N0(G,  £)  where  iV0  is  linear.  If  one  replaces  G  =  (y+(x,y), ....) 
from  (3.6)  by  G(ft,p)  =  (yi(V(Qi(x)t(y)  -f  i?(x,  y))+, ...)  from  (4.3)-(4.4)  and  defines 
Ar(ft,£,p)  to  be  iVo(G(ft,p),  £),  one  obtains  a  nonlinear  equation  ft  =  JV(ft,£,p)  whose 
fixed  points  formally  provide  solutions  of  (4.1)-(4.2). 

As  is  standard  in  an  analysis  leading  small  solutions  in  a  center  manifold  we  shall 
use  a  truncation.  Let  77  :  R  — ►  R  be  an  odd,  nondecreasing,  G°°  function  which  is  the 
identity  for  —  1  <  t  <  1,  equal  to  2  for  t  >  2,  and  satisfies  y'  <  2.  Set  r/r(t)  =  rr](t/r )  for 
r  >  0.  The  terms  g\,  g2,  and  g  in  (4.3)-(4.4)  can  be  expressed  in  terms  of  Q,  Q',  R  and 
their  derivatives,  or  in  terms  of  Q 1,  Q 2,  R,  and  their  derivatives.  In  such  an  expression 
we  replace  each  occurence  of  these  latter  variables  by  its  composition  with  rjr.  Thus,  for 
example,  g 1  =  wxwy/{\  +  wy)  is  replaced  by 


(4.6) 


9  i(ft,ri  0) 


{pT]r(Q2(x))t(y)  +  T?r(Rr(^,!/))}{^r((?l(x))f/(y)  +  T)r(Ry(x,  y))} 
1  +  Vr(Ql(x))t'(y)  +  rjr(Ry(x ,  y)) 


and  similar  expressions  g2(Cl,r,l3),  g(Q.,p,r,fi)  are  defined.  Now,  with  G  in  equation  (4.5) 
replaced  by  G(ft)  =  G(ft,p,  r,/3)  =  (gf(n,r,0),...,g(U,p,r,0))  one  obtains  an  equation 

(4.7)  ft  =  iV(ft,  p,  r,  /?) 
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where  N  has  components 

Ni(to,$,P,r,0)  =6  +  /?  f  rfr(Q2)(s)ds , 

Jo 

(4'8)  Nt(Q,(,P,r,ff)  =  (t  +  ff-lLlGlQ) IJ  +  /T1  [‘ L2G(Sl)(s)ds, 

Jo 

N3(n,t,P,r,P)  =  ( L3G)(x,y ),  (x,y)  G  T. 

Before  proceeding  to  an  existence  result  for  (4.7),  some  remarks  are  in  order.  A 
solution  u>  corresponding  to  a  wave  will  be  small  in  Xq  ;  that  is,  w  and  its  derivatives  will 
be  uniformly  small.  However,  the  operators  N\  and  N2  do  not  map  bounded  functions  into 
bounded  functions.  They  do  map  spaces  with  growth  parameter  p  >  0  into  themselves, 
and  that  is  where  we  are  obliged  to  begin.  In  [3]  operators  similar  to  the  Ni  arose  through 
a  solution  process  involving  Fourier  transforms.  There,  however,  the  maps  had  sufficient 
regularizing  properties  that  it  was  not  necessary  to  introduce  the  ’small’  parameter  ft  nor 
to  introduce  the  sets  used  here.  In  the  ensuing  subsections  we  shall  repeatedly  use 
results  from  [6]  regarding  the  behavior  of  nonlinear  maps  in  weighted  spaces. 

4.2.  Existence  with  truncation 

The  basic  nonlinear  map  of  interest  is  the  composition  of  an  element  u  =  (u1? ...,  un) 
with  a  smooth  function  g  from  Rn  to  R.  More  generally  we  consider  a  smooth  function  g  = 
g(ui,  ...,un,pi,  depending  on  a  parameter  p  —  (pi,...,pn),  though  the  dependence 

on  the  parameter  p  is  not  central  in  the  discussion  and  will  sometimes  be  suppressed.  We 
use  Dk  for  a  kth  derivative  of  a  map  between  Banach  spaces  and  |Pfc(-)|  for  its  norm 
(cf.  [1]).  It  will  be  assumed  that  the  derivatives  with  respect  to  u  satisfy  < 

oo,  k  =  0, 1,2, ..  for  all  (u,p)  €  Rn  x  Rn.  We  use  C  for  a  generic  constant  depending  on 
j,  n,  and  other  quantities  indicated  in  the  context. 

Suppose  g  is  smooth  and  that  y(0,  ..0,pi,  ..,p„)  =  0.  Let  Mg  be  the  map  defined  by 

(4.9)  Mg(u,p)(x,y)  =  g(u(x,y),p). 

If  p  >  0,  then  Mg  maps  C^“(T, R")  x  R"  into  C^’%  where  a  <  Csup l<t<J+lckSk.  Of 
course,  we  shall  be  interested  in  truncated  functions  as,  for  example,  in  (4.6).  We  define 
y(u,p,r)  to  be  y(pr(wi), ...,  pr(wn),Pi,  ■~,Pn)-  For  such  a  composition  g,  results  regarding 
Mg  needed  for  an  existence  theorem  are  now  summarized  (cf.  [6]). 
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Proposition  4.1.  Suppose  h(u,p )  is  smooth  and  that  h(0,p )  =  0.  Let  g  denote  the 
truncated  function  h{u,p,r).  Then  Mg  maps  Rn)  x  Rn  into  with 

<7  <RCr(l  +  |p|), 

where  B  =  SUP  1<k<j+1  ($/r)k  and  C  depends  on  h.  Further,  \Mg\ioo  <  Cr  and 

(4.10)  | Mg(u,p)  -  Mg(u',p')\  <  BC(\u  -  u'\  +  r(p  -  p' |), 

where  the  norms  of  functions  are  those  in  Cf.  If,  in  addition,  g  is  flat  (i.e.,  satisfies 
Du,Plg{ 0,  0, 0,p2,  pn )  =  0  for  all  p),  then  Mg  maps  Rn)  x  R"  into  Cfc*  where 

a  <  BCr(r  +  |pi|). 

Further,  <  Cr(r  +  |pi|)  and 

(4.11)  | Mg(u,p)  -  Mg(u',p')\  <  BC(r  +  |pi|)|u  -  u'l  +  RCr|p  -p'|. 

Because  of  the  explicit  nature  of  the  map  Ni  we  can,  and  will,  choose  6  =  r  in  what 
follows,  so  that  B  ~  1.  Recall  that  lies  in  17 ^  =  Ci’,r^2(R)  x  C^’,r^2(R)  x  X*  r.  In  (4.8), 
p  =  pi  =  Xd  —  X  and  we  drop  the  subscript.  To  see  that  N  maps  into  itself  for  small  r 
and  is  contractive,  we  can  ignore  the  additive  terms  £i, £2-  The  remaining  terms  all  have 
the  form  of  compositions:  linear  maps  (derivatives  and  injections),  followed  by  nonlinear 
functions  of  the  type  covered  by  Proposition  4.1,  followed  by  other  linear  maps.  We  use 
pairs  or  triples  of  numbers  to  indicate  the  sup  remum  bounds  on  a  function.  Thus  for  N 1  in 
(4.8),  from  Proposition  4.1  or  by  direct  computation,  it  follows  that  the  map  Q 2  — >  Pr(<?2) 
takes  a  function  with  bounds  (00,  r)  on  the  function  and  its  Holder  estimate,  respectively, 
to  one  with  bounds  (2  r,  2r).  The  linear  map  /3  f*  produces  bounds  (00, 2  /3r,  2  fir),  the  third 
term  indicating  the  Holder  bound  on  the  first  derivative.  The  Lipschitz  estimate  is 

1^(0)  -  <  t3C\Sl  - 

where,  because  of  the  integration,  C  depends  inversely  on  p. 

In  iV2  the  linear  maps  d/dx  and  d/dy  take  a  function  with  bounds  (00,  r,  r)  to  one  with 
bounds  (r,  r).  The  injection  of  C]fX^2  into  C^’1^2  is  a  bounded  linear  map  and  produces 
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functions  with  bounds  (oo,r).  Hence  the  nonlinear  flat  functions  in  G(ft)  act  on  functions 
with  bounds  (oo,r)  and  produce  (Cr(r  +  |p|),(Cr(r  +  |p|),  according  to  Proposition  4.1. 
Both  L\  and  L2  preserve  such  bounds  with  a  new  constant  C.  The  integration  following 
£2  gives  (00,  Cr{r  +  |p|)).  In  all,  iV2  produces  bounds  (00,  /3~1Cr(r  -f  |p|))  and,  from  (4.11), 
satisfies 

(4.12)  |JV„(n,p)  -  <  l[C(r  +  |p|)|fi  -  n'lc,  +  Crlp-p'll 

In  JV3,  as  earlier,  the  term  G(Q)  produces  terms  with  bounds  (Cr(r  +  |p|),(Cr(r  +  |p|)). 
As  L 3  acts  from  A0  to  Xq  (cf.  Remark  3.4),  there  is  a  bound  C'r(r  +  |p|)  for  the  resulting 
function,  its  derivatives,  and  their  Holder  estimates.  The  right-hand  side  of  the  Lipschitz 
estimate  for  in  X *  is  the  same  as  that  for  jV2,  but  without  the  factor  /?-1  so  that 

(4.13)  |jV(Sl,p)-JV(SlV)k  <  C,(/J  +  (l+/J-,)(r+|p|))|n-n'|t,.+C,(l  +  )J-1)r|l>-p'l- 


In  order  to  carry  out  a  contractive  scheme  for  N  we  require  (i)  r  6  (0,  ry]  so  that  all 
denominators  in  terms  gi(u,p,r),  etc.  are  at  least  1/2;  (ii)  2/?  <  1  and  Ci/?-1(r  +  |p|)  <  1 
so  that  N  maps  into  U^y  and  (iii)  C2(/?  +  (l  +  /3~1)(r  +  |p|))  <  1  so  that  it  is  contractive  in 
the  argument  ft.  Here  C\  and  C2  arise  from  combining  terms  from  the  last  paragraph.  To 
satisfy  these  conditions  with  a  contraction  factor  of  1/2,  first  choose  /?  =  min((4C2)-1, 1/2) 
and  restrict  p  by  jp|  <  r.  Then  there  exists  an  vq  >  0  such  that  (i)-(iii)  hold  for  any 
r  6  (0,  r 0] - 

Theorem  4.2.  There  exist  a  positive  constants  f3  and  r0  such  that  if  r  6  (0,ro],  the 
equation  Q  =  N(S2,  £,p,  r,  jd)  has  a  unique  solution  Q  G  I7^,r  for  each  choice  of  £  E  R2  and 
p  =  Xd  —  A  G  (— r,  r).  Moreover,  is  Lipschitz  continuous  as  a  function  of  £  and  p. 

We  fix  r  =  r0  and  (3  from  the  theorem  and  suppress  them  when  convenient.  Let 
ft(£,p)  =  ft(£iP)(z,y)  =  (Qi(Z,p)(x),Q2(Z,p)(x),R(Z,p)(x,y))  denote  the  solution  given 
by  Theorem  4.2. 

Remark  4.3.  If  £1  =  £2  =  0,  then  the  trivial  solution  satisfies  the  fixed  point  equation, 
whence  ft(0,p)  =  0,  p  G  (— r,  r). 

Remark  4.4.  The  constants  C\  and  C2  in  the  discussion  preceding  Theorem  4.2  depend 
on  p  and  will  approach  00  as  p  — >  0.  Suppose  I  =  [pi,p2]  where  0  <  pi  <  p2  <  p  and  let 
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Ci,  i  =  1,2  be  the  supremum  of  Ci{p)  for  p  G  I.  A  choice  of  (3  and  ro  based  on  these  new 
constants  provides  a  solution  for  any  p  £  I  and,  by  uniqueness,  a  solution  in  U ^  coincides 
with  that  for  any  larger  parameter  in  /. 

4.3.  Reduction  and  regularity 

Here  the  special  structure  of  the  third  component,  R(x,y ),  is  examined.  It  is  the 
dependence  on  x  that  is  central  here  and  we  can  naturally  suppress  the  variable  y  by 
regarding  R(x,  •),  Rx(x,-),  and  Ry(x,-)  as  talcing  values,  respectively,  in 

v0  =  {feC(i): 

Vi  =  {/  €  C(I)  :  f  €  C0-1/2^)}, 

v2  =  {/  e  c,0’1/2(7±)}. 

Let  xo  be  fixed.  As  for  ordinary  differential  equations,  one  defines  fl*(x)  =  fi(£,p)(x  + 
xo),  fi**(x)  =  fl(Qi(£,p)(x0), Q2(£,p)(xo),p)(z)  and  shows  (cf.  [3], [6])  that  they  satisfy 

fi*(x)  -  0**(x)  =  N(Q*(x))  -  N(Q**(x)). 

As  they  have  the  same  data  at  x  =  0,  they  coincide,  yielding  the  following  important 
result. 

Lemma  4.5.  For  a 11  £  €  R2,  p  €  (— r,  r),  x  €  R,  and  xo  E  R 

ft(£»P)(z  +  x0)  =  S2(Qi(£,p)(z0),  Q2(C,p)(*o),p)(z). 

By  taking  derivatives  with  respect  to  x  and  y  in  the  third  component  of  this  last 
relation,  setting  x  =  0,  and  then  replacing  xo  by  x  one  obtains  the  local  dependence  of  R 
and  its  derivatives  on  Q\{x)  and  Q2(x)  expressed  in  the  following  ’center  manifold’  result. 

Theorem  4.6.  There  are  Lipschitz  continuous  functions  J0,  J\,  and  J2  from  R2  x  (— r, r) 
into  Vo,  V\,  and  V2,  respectively,  such  that 

ft($,p)(x,-)  =  MQi(CpXx),Q2(^,p)(x),p) 

(4-14)  Rx($,p)(x,  •)  =  Ji(Qi((,p)(x),  Q2(Z,p)(x),p) 

Ry((,p)(xw)  =  J2(Qi(Z,p)(x),  Q2(Z,p)(x),p) 

At  this  point  we  indicate  a  difficulty  that  arises  in  the  context  of  weighted  spaces, 
where  the  analysis  eventually  requires  some  degree  of  smoothness  in  the  resulting  solutions. 
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Let  g(u )  =  T]l(u )  for  r  >  0.  As  a  map  of  C°  into  itself  the  map  Mg  is  not  differentiable  at 
u  =  0.  Were  there  a  derivative,  it  would  take  each  uniformly  bounded  function  in  C°  to 
zero.  However,  if  h  £  C°(R)  is  a  function  with  compact  support,  the  function  h *  =  h(x  —  k ) 
has  C\  norm  \hk\  <  Ce~k  — *  0  as  k  — >  oo.  However,  one  verifies  that  the  norm  \ql(hk)\ 
does  not  approach  zero  faster  than  \hk\.  One  can  show  that  Mg  is  Lipschitz  continuous 
from  Ci  into  itself  (cf.  Remark  4.8  below).  However,  it  is  not  Lipschitz  from  C]'a  into 
itself  for  p  >  0  when  either  j  or  a  is  positive.  The  lack  of  regularity  of  Mg  is  attributable 
to  the  presence  of  multiplicative  factors  which  can  have  exponential  growth  in  expressions 
for  derivatives.  One  can  control  these  factors  if  one  views  the  target  of  Mg  as  a  space 
allowing  more  rapid  growth  than  that  of  the  domain  and  for  applications,  differentiability 
in  this  sense  will  suffice. 

For  fixed  (x,  y)  €  T  and  q  =  1,2,..  Dqg  =  D^u  p.)g(u(x,y),p )  is  a  symmetric  q-linear 
map  from  R”  x  R”  to  R.  It  induces  a  symmetric  q-linear  map  on  Ci,a  x  Rn  defined  by 

MD«g{hi ,  7TJ )....{hq,nq)  =  Dqg(u(x,y),p)(h1(x,  y ),  ttj )..{hq{x,  y),  irq). 


The  following  result  is  shown  in  [6]. 

Theorem  4.7.  Suppose  g  is  smooth  and  p  >  0.  Then  Mg  defined  in  (4.9)  is  q  times 
differentiable  as  a  map  from  C^Q{T,  Rn)  X  Rn  to  C^-“  ?+2)M  372 the  derivatives,  defined 
by 

Dk  Mg  =  MDkg ,  1  <k<q 

are  Lipschitz  continuous. 


Remark  4.8.  While  DMg  does  not  generally  exist  as  a  map  of  C2’"  into  itself,  the  map 
Mog  is  a  bounded  linear  map  of  C2,a  into  itself  when  evaluated  at  a  u  G  Cp° 

From  Theorem  4.7  one  can  obtain: 

Corollary  4.9.  For  p  >  0  the  map  N(Q,  £,  p)  has  q  Lipschitz  continuous  derivatives  acting 
from  Up  x  R2  x  R  to  U(q+2)n- 

In  the  spirit  of  Remark  4.8  one  can  show  that  the  formal  derivative  of  N  with  respect 
to  fl,  evaluated  at  a  solution,  is  a  bounded  map  of  U ^  into  itself  with  a  bound  equal  to  the 
Lipschitz  constant  in  (4.13),  at  most  1/2  under  our  working  hypotheses.  It  is  known  [11] 
that  functions  describing  center  manifolds,  such  as  those  in  Theorem  4.6,  are  not  C°°  in 
general.  However,  using  Remark  4.4  and  the  previous  corollary  it  can  be  show  inductively 
(cf.  [6])  starting  with  a  sufficiently  small  p  that  the  following  holds. 
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Theorem  4.10.  Given  q  >  1  there  is  an  integer  k(q)  such  that  if  p  =  k(q)p  <  ~p,  then 
the  solution  Q  €  UM>r  from  Theorem  4.2  is  q  times  differentiable  with  respect  to  (£,  p)  as 
a  map  into  Up. 

A  final  result  will  be  needed  for  the  analysis  in  section  5. 

Theorem  4.11.  Under  the  hypotheses  of  the  previous  theorem  the  functions  J,,  i  =  0, 1, 2 
in  Theorem  4.6  are  q  times  differentiable  and  satisfy 

(4.15)  \Ji^p)\v,  <CM|£[  +  |p(). 


Proof.  Let  Eq  be  the  evaluation  map  taking  R,  Rx,  or  Ry  in  X etc.  to  its  value  at 
x  —  0  in  Vo,  Vj,  or  V2i  respectively.  Since  Eq  is  linear  and  hence  infinitely  differentiable, 
the  smoothness  of  the  J{  follows  from  Theorem  4.10. 

One  has  D^_pfl  =  (7  —  D^N)~l  D^tPN  and  calculates  that  D&N  evaluated  at  Cl  — 
0,  p  =  0  takes  (ux,U2,u3)  to  (/?  f*  U2(s)ds,  0, 0).  Since  neither  fi  nor  £ 2  appears  in  N3 
and  since  p  intervenes  through  flat  functions,  one  finds 

D{li^,pft(0,0,0)(a1,or2,7r)  =  (CTi  +  ^<r2x,  ct2, 0) 

and  thus  the  third  component  is  the  zero  function.  It  follows  that  the  derivatives  of  i?, 
Rz ,  and  Ry  with  respect  to  £,p  vanish  at  (0,0).  From  Remark  4.3  each  of  these  functions 
vanishes  identically  when  £  =  0  and  p  €  (— r,  r)  so  that,  for  example,  R(£,p)  =  S(£,p)£. 
Since  Z)^R(0,0)  =  0,  5(0,0)  =  0,  implying  that  5  vanishes  at  least  linearly  in  |£|  +  \p\. 
Similar  considerations  apply  to  Rx  and  Ry,  yielding  the  desired  bounds  on  Jj,  i  =  0, 1,2. 
q.e.d. 

5.  THE  REDUCED  EQUATIONS 
5.1.  Structure  of  the  fixed-points 

In  this  section  we  analyze  solutions  of  (4.7)  for  A  in  a  neighborhood  of  Aj  and  it  will 
be  convenient  to  express  all  quantities  in  terms  of  p  =  A<f  —  A  throughout  most  of  the 
section.  For  fixed  £  and  p,  write  fi(x,y)  for  fi(£i,£2 ,p){x,y),  f2*(x,y)  for  fl(£i ,  —  £2)(z>  y) 
and  let 


n(i,y)  =  (Qi(£i,  -£2,p)(-*),  -Q2(ti,-h,p)(-x),R(£i, 


-$2,p)(-*,y))- 
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The  particular  form  of  the  truncated  g\  in  (4.6)  and  its  analogues  for  g 2  and  g  give  the 
relations  gi(&)(x,y)  =  -g i(ft*)(-x,y),  g2(Cl)(x,y)  =  g2(Ct*)(-x,y),  and  g(Q)(x,y)  = 
y(fl*)(— x,  y).  Use  of  the  explicit  forms  of  and  L2  from  (3.30)-(3.31)  in  (4.8)  yields 
N t(Q. ,  ,  (2 ,  p)(x)  =  f2,(x)  for  i  =  1,2,  while  the  linear  analysis  of  Section  3  (cf.  (3.6)- 

(3.7))  proves  the  analogous  relation  for  iV3.  Hence,  N(Q,  £1, &,p)(x)  =  ft(x).  Since  D 
satisfies  the  same  equation  and  N  is  a  contraction,  there  results  fl(x,y)  =  fi(x,y),  giving 
the  following  useful  structure  relations: 

Qi(6,6,P)(*)  =  Qi(Zi,-Z2,p)(-x), 

Q2ttl,t2,P)(x)  =  -Q2ttl,-t2:P)(-x), 

(5.1)  R(ti,t2,p)(x,y)  =  R(Zi,-b,p)(-x,y), 

RxUi,(2,p)(x,y)  =  -Rx(Zu-b,p)(-x,y), 

Ry(Zl,b,P)(x,y)  =  Ry(tu-Z2,p)(-X,y)* 

for  all  (^1,^2)  €  R2,  P  =  \d  -  A  6  (-ro,r0),  and  (x,y)  e  T. 

5.2.  A  differential  equation  for  Q 

Let  s,  s  1  €  (0,ro),  where  r0  arises  in  Theorem  4.2,  and  assume  that  for  some  (6R2 
and  p  €  [— s,s]  there  holds 

(5-2)  IQi(^p)(OIlco(R-)>  \Q2(Z,p)(-)\l00(K)  <  si- 

The  use  of  this  in  (4.15)  gives  the  following  estimate  whenever  Si  is  sufficiently  small: 

(5-3)  l-R(£;p)(-)U«,mi  I#i(£;p)(-)Uco(t)>  l-Ry(f;p)(-)Uoo(T)  <  r0. 

These  estimates  allow  one  to  remove  the  truncations  present  in  (4.8).  Write  Q(x)  for 
Qi(^;A)(x),  so  that  (3~lQ'(x)  =  Q2(t,p)(x)  by  (4.5).  Differentiating  the  second  equation 
in  (4.8)  yields 

(5.4)  Q"(x)  =  pj-iQ2(f,,p){*)  =  ^I,(G(f!(Cp)))(*)  +  £2(C(S!«,p)))(x). 

Wherever  R,RZ,  or  Rv  occurs  in  the  right-hand  side,  we  use  (4.14)  to  replace  it  by  a 
function  of  Q(x)  and  Q'{x).  There  results 

(5.5)  Q"(x)  =  (Q(x),Q'(x),p)  +  E2(Q(x),Q'(x),p) 

ax 
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where  E\  and  E2  are  smooth  functions  of  their  arguments,  by  Theorem  4.11.  One  may 
check  easily  that  Ei  and  its  gradient  vanishes  at  (0,0,0)  for  i  =  1,2.  If  s  and  si  are 
sufficiently  small,  then  by  solving  for  Q"  in  the  equations  resulting  from  (5.5)  we  are  led 
to  an  equation  of  the  form 

Q"(x)  =  F(Q(x),Q'(x),p),  xeR 

(5.6) 

Q(0)  =  6,  Q'( 0)  =  ^2- 

We  have  shown  that  if  (5.2)  holds  for  some  suitably  small  si,  then  a  solution  Q(x)  = 
Qi(£,p)(x)  to  (5.6)  exists.  Conversely,  if  one  can  construct  a  solution  (p,  Q)  of  (5.6)  with 

(5.7)  |<3Il00(R),  0~1\Q'\loo(K)  <  -si,  bl  <  5, 

then  setting  Qi((i,p)(x )  =  Q(x),  Q2(£,p)(x)  =  P^Q'ix),  and  (cf.  (4.14)) 

(5.8)  R(x,y)  =  R(£,p)(x,  y)  =  R(Q(x),  Q'(x),p)(0,  y) 

yields  Q (£,p)  satisfying  (4.8)  without  any  truncation.  In  particular,  setting  w(x,y)  = 
Q(x)t(y)  +  R(x,y )  yields  a  small  solution  of  (1.11)-(1.13).  Hence,  the  study  of  small 
solutions  to  (1.11 )-( 1  - 13)  for  p  =  \d  -  A  small  is  reduced  to  studying  (5.6). 

The  structure  results  (5.1)  easily  show  that  Ei  is  an  odd  function  of  its  second  argu¬ 
ment  Q'(x )  while  E2  is  even.  It  follows  immediately  that  F  in  (5.6)  is  an  even  function  of 
its  argument  Q'(x).  This  ensures  that  any  solution  of  (5.6)  is  necessarily  symmetric  about 
a  point  at  which  Q'  vanishes,  whence 

Theorem  5.1.  Let  Q  satisfy  (5.6)-(5.7).  Then  exactly  one  of  the  following  holds:  (i)  Q{x) 
=  constant  ( Q'(x )  =  0);  (ii)  Q  is  monotone  (Q'(x)  ^  0  on  R);  (Hi)  Q  is  even  about  some 
point  xo  and  is  monotone  on  (r0,  oc)  (Q'(x)  has  a  single  zero  on  R);  (iv)  Q  is  a  periodic 
function  (Q1  has  more  than  one  zero  on  R). 

This  result  has  implications  for  any  small  solution  w.  Since  w(x,y)  =  Q(x)t(y )  + 
R{x,y),  where  R(x,y)  =  R(Q(x),  /?-1Q'(r),p)(0,  y),  there  results  wz{x,y)  =  Q'(x)t(y)  + 
Rz(Q(x),l3~lQ'(x),p)( 0,y)  by  (4.14).  This  may  be  rewritten  as 

(5.9)  u;I(x,y)  =  Q'(x){t(y)  +  Q'(^)"1^(Q(^)>/?_1^V),p)(0,y)}. 

Now  Rx(£i, £2;p)(0,  y)  is  an  odd  function  of  £2  by  (5.1),  and  so  ,  f,2,p)(®,  y)  is 

an  even  function  of  £2  and  is  smooth  by  Theorem  4.11.  This  quantity  vanishes  in  Vj  when 
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£i  =  £2  =  0  and  p  =  0.  It  follows  that  the  expression  added  to  t(y )  in  the  brackets  in  (5.9) 
is  small  in  Vj,  uniformly  in  x,  whenever  s  and  si  are  sufficiently  small.  Elliptic  theory 
then  shows  the  added  term  is  small  in  a  topology  including  derivatives  near  y  =  —  h  and 
y  =  1  —  h.  Thus  the  whole  bracketed  term  in  (5.9)  is  positive  on  T  whenever  s  and  Si 
are  sufficiently  small.  Theorem  5.1  then  ensures  that  w  satisfies  one  of  the  following:  (i) 
w(x,y)  is  a  function  of  y  alone;  (ii)  wz(x,y )  is  one-signed  on  T;  (iii)  u>(x,y )  is  even  in 
x  about  some  x0  and  wx(x,y)  is  one-signed  on  (x0,oo)  x  7;  (iv)  to  is  a  periodic  function 
of  x.  In  case  (i),  w  is  either  identically  zero  or  is  a  nontrivial  x-independent  solution  of 
((1.1 1)-(  1.13);  these  are  termed  conjugate  solutions  by  Benjamin  [7], [8].  In  case  (ii),  we 
have  a  heteroclinic  orbit  connecting  a  conjugate  flow  at  x  =  — oo  to  a  different  one  at 
x  =  +oo.  We  shall  show  there  are  solutions  connecting  the  trivial  solution  to  a  non-trivial 
conjugate  flow;  these  are  ‘internal  surges’.  After  a  suitable  translation,  we  may  take  xo  =  0 
for  case  (iii);  then  the  flow  connects  a  conjugate  flow  at  x  =  — oo  to  itself  at  x  =  -f  oo,  and 
so  is  a  homoclinic  orbit.  If  the  conjugate  flow  is  the  trivial  solution,  then  the  solutions  are 
usually  called  solitary  waves. 

One  may  seek  small  solutions  of  (5.6)-(5.7)  by  expanding  F  to  quadratic  or  cubic 
terms  in  ( Q,Q',p ),  scaling  the  resulting  equations,  and  analyzing  them  with  the  aid  of  an 
implicit  function  theorem.  This  approach  is  used  in  [3]  for  a  water-wave  problem  and  in 
[5]  for  the  case  of  internal  waves.  One  may  avoid  these  technicalities  with  the  use  of  an 
important  conserved  quantity,  and  this  is  now  shown. 

5.3.  A  first  integral 

As  noted  in  (1.17)  or  [7],  the  flow  force  is  independent  of  x.  In  terms  of  the  function 
w,  this  becomes 

r 

(5.10)  /  Poo-^ - z-dy  +  A (p  -  l)u?2(x,  0)  =  C 

Jr  (l+Wy) 

for  all  x,  where  C  is  a  constant.  This  result  can  also  be  derived  from  (1.11 )-( 1.13): 
one  multiplies  (1.11)  by  poowz  arid  integrates  over  (  —  h,  1  —  h )  using  (1.12)-(1.13).  Since 
wy(x,y)  =  Q(x)t'(y)  -1-  Ry(Q(x), (3~1Q'(x),p)( 0,y),  it  follows  from  (5.1)  that  wy  is  an 
even  function  in  the  argument  Q'(x).  By  Theorem  4.6,  we  may  assume  that  u?(x,y)  and 
wy{x,y )  are  functions  of  Q(x),Q'(x)2  and  y.  Equation  (5.9)  gives  tx;2(x,y)  as  a  function 
of  Q(x),  Q'(x)2,  and  y,  which  behaves  like  Q'(x)2t2(y)  if  s  and  Sj  are  sufficiently  small. 
Hence,  (5.10)  has  the  form 

(5.11)  H(Q'(x)2  ,Q(x),p)  =  C 
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where  dH/dai(ai,a.2,p)  has  the  value: 

Poo(y)t2(y)dy  =  -1 

when  =  <32  =  0  and  p  —  0.  To  prove  this  one  uses  (4.15)  in  (5.10).  It  follows  from  (5.11) 
with  the  aid  of  the  implicit  function  theorem  that 

(5.12)  Q'{x)2  =  h(Q(x),p),  xeR, 


where  h  depends  on  the  constant  C.  Recall  that  F  in  (5.6)  is  an  even  function  of  Q'\  as 
F  is  smooth,  it  is  a  function  of  ( Q ')2.  The  use  of  (5.12)  in  (5.6)  yields 

(5.13)  Q"(x)  =  F(Q(x),p). 


The  function  k  in  (5.12)  is  related  to  the  primitive  of  F. 


5.4.  Small  solutions  of  equations  (5.12)-(5.13) 

One  may  seek  solutions  of  (5.12)-(5.13)  satisfying  (5.7),  depending  on  the  parameters 
p  and  C.  Since  the  original  interest  of  this  paper  was  with  internal  surges  (solutions 
connecting  the  trivial  solution  at  x  =  — oo  to  a  non-trivial  conjugate  flow  at  x  =  +oo),  we 
begin  with  the  case  Q(x)  — *  0  as  x  — ►  — oo,  whence  C  =  0  in  (5.10).  Periodic  solutions  will 
be  found  later  by  a  phase-plane  argument. 

Before  the  differential  equation  for  Q  is  studied,  a  few  remarks  are  necessary  about 
the  parameters  in  the  problem.  The  numbers  p  and  h  have  been  fixed  in  the  interval  (0, 1). 
All  estimates,  such  as  those  needed  for  the  contraction  mapping  hold  uniformly  for  p  and 
h  varying  in  compact  subsets  of  (0, 1).  Instead  of  working  with  p  and  h,  we  take  e  and  p 
as  independent  parameters  and  define  h(e,p)  implicitly  through  (cf.  (1.15)) 


(5.14) 


_1 _ P__ 

6  h 2  (1  -hf 


When  convenient,  we  shall  continue  to  write  h  rather  than  the  function  of  e  and  p. 
Given  a  fixed  positive  number  k ,  we  shall  restrict  attention  to  e  6  [—&,&]  and  p  with 
p~l  <  k ,  jl  —  p\~l  <  k  .  The  use  of  the  representation  w(x,y )  =  Q(x)t(y)  +  R(x,  y),  where 
R(x,y)  =  R(Q(x),  fl~1Q'(x),p)(0,y),  together  with  (4.15)  in  (5.10)  with  C  =  0  yields 


(5.15)  Q'(x)2  =  t(0)2p(l  -p)Q2(x)  +  0(pQ2(x))  +  0(Q3(x)),  x6  R 
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as  Q,p  — >  0.  As  noted  earlier,  i?(£i,£2,p)(0,y)  is  a  smooth  function  of  (£i,£2,p)  which  is 
even  in  £ 2 •  In  addition,  i?(0, 0; p)(0,  y)  =  0  for  all  p  6  (— ro,ro)  and  the  Frechet  derivative 
of  the  map  (£i,£2,p)  — ►  i?(£i,  £2,p)(0,  y)  vanishes  at  £1  =  £2  =  P  =  0.  The  use  of  this 
together  with  (5.15)  gives 

(5.16)  R(x,y)  =  pa(y,p)Q(x )  +  b(y,p)Q2(x)  +  Q3(x)d(Q(x),y,p), 


where  a,  b,  and  d  are  smooth  functions  of  their  arguments.  One  can  use  this  representation 
in  (5.4),  and  after  lengthy,  but  straightforward  calculations,  arrive  at  the  equation 


(5,17) 


3 

Q"(x )  =  pai(p,  e)Q(x)  +  -(ea2(p,  e)  +  pa3(p,  e))Q2(x) 
+  2a4(p,  e)Q3(x)  +  Q4(x)a5(Q(x),p,  e) 


where  ai, ...,  as  are  smooth  functions  of  their  arguments.  These  functions  also  depend  on 
p,  but  the  dependence  is  suppressed.  The  following  holds  for  the  coefficients  in  (5.17): 


(5.18) 


ai(0,e)  =  m2h2(  1  —  p), 
a2(0,  e)  =  —m3h3, 

«3(0,0)  =  ^~~~m5h5 

a4(0,  0 )  —  m  h  (  —  +  ^  _  ^3) 


where  m  is  the  normalization  constant  in  (3.9): 

2  _ 3 _ 

777  —  - 

h*(h  +  p(i-h)y 

Our  task  is  to  find  the  totality  of  solutions  (p,  Q)  to  (5.17)  as  e  varies  over  [—it,  k]  with 


(5.19) 


Q(x)  — »  0  as  x  — >  —00, 


subject  to  the  constraints  (cf.  (5.7)) 

(5.20)  |p|  <  s  and  0  <  |<5|Z/oo(R)  <  sj. 


The  condition  in  (5.7)  that  |(3'(:r)|  <  is  implied  by  (5.15)  and  (5.20)  whenever  s 
and  Sj  are  sufficiently  small. 
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If  (5.17)  is  multiplied  by  Q'  and  integrated  with  the  aid  of  (5.19)  then 
(5.21)  e'(i)2  =  02Wn(Q(x),p,e), 

where 


II(<2,p,  e)  =  pax(p,  e )  +  (ea2(p,  e)  +  pa3(p,  e))Q  +  a4(p,  e)Q2  +  Q3a5(Q,  p,  e) 


and  5,5  is  related  to  a5  from  (5.17)  in  an  obvious  way. 

The  next  step  is  to  find  the  zeros  of  II (•,  p,  e)  and  their  dependence  on  (p,  e).  Note 

that 

11(0, 0,  e)  =  0,  e  €  [—  k,  Ar], 
np(0, 0,  e)  =  aifO,  e)  >  0,  e  £  [— Ar,  A:], 

(5-22) 

nQ(0, 0,  e)  =  ea2(0,  e )  =  —m3h3e  ^  0  for  e  7^  0, 

nQQ(0,0,0)  =  2a4(0,0)>0. 


Assume  for  the  moment  that  there  is  a  triple  (Q,p,  e)  €  [— Si,Si]  x  [— s,s]  x  [— k,k] 
for  which  II  =  IIq  =  0.  Then  Q(-)  =  Q  satisfies  (5.17)  and  (5.21),  whence  w(x,y)  = 
Qt(y)  +  R(Q,Q,p)  is  a  solution  of  (1.11)-(1.13)  satisfying  (5.10)  with  C  =  0.  Such  solutions 
were  considered  in  section  8  of  [4]  and  the  following  shown: 

A_(l  +  yp)2 

1  -  p 

(5.23)  p (e>  p)=*d-*= X(T^7)(1 "  ^i~eh2y 

qu  p )  =  — _ £ _ s  vgygyrrgg 

m(l  +  yp)(i  +  ^)  6  (l  +  y'p)9/2 


where  h  is  considered  a  function  of  p  and  e  (cf.  (5.14)).  Since  |p|  =  |A^  —  A|  <  s,  and  s  is 
to  be  small,  we  have  that  e  is  small  and 


(5.24) 


P(e,p)  =  [7 


1  y/p 


4(1  p)  ( 1  +  s/p\ 


r]e; 


so  that  Q  is  linear  in  e  and  p  is  quadratic,  to  lowest  order. 

Since  the  function  II(Q,p,  e)  has  IIp(0,  0,  e)  >  0,  by  (5.18)  and  (5.22),  there  is  a  surface 
S  given  by  a  function  p  =  P{Q,  e)  for  ( Q ,  e)  in  a  neighborhood  of  (0, 0),  which  is  precisely 
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the  locus  of  zeros  of  II  in  a  neighborhood  of  (0,0,0).  The  surface  contains  the  curve  of 
conjugate  solutions 

r  =  {(Q(e),p{e),e)  :  —k'  <  e  <  k'} 

for  some  0  <  k'  <  k,  where  now  we  suppress  the  dependence  on  p.  On  T  one  has  II  =  0, 
IIq  =  0,  and  IIqq  0  from  (5.22).  It  follows  that  in  a  neighborhood  of  (0,0,0)  the 
following  derivative  relations  hold:  Pq  =  0  on  T  and  Pqq  =  -Uqq/Up  <  0  on  S.  Hence 
the  surface  5,  viewed  in  terms  of  p  and  e  as  parameters,  is  composed  of  two  smooth  sheets, 
Q  =  Qi(Pi e )  lying  above  Q  =  Q2(p,  e ),  and  the  curve  T  along  which  they  meet  and  become 
singular.  Since  Q(e)  has  the  same  sign  as  e,  the  upper  sheet  is  distinguished  for  p  =  0 
by  having  Qi  =  0  for  e  <  0,  while  Q2  =  0  for  p  =  0  and  e  >  0.  Since  Hq(0,  0,  e)  ^  0 
for  e  7^  0,  it  is  also  clear  that  the  distinguished  sheet  can  be  continued  to  nonzero  p  in  a. 
neighborhood  of  each  e  ^  0.  To  summarize  these  results  in  a.  convenient  form  we  introduce 
the  following  notation:  For  sufficiently  small  ,s,  let  ei  €  ( —k,  0),  e2  €  ( 0,  A- )  Ire  the  small 
roots  of  s  =  p(  • ,  p ) ,  so  that 

C<  y/P 

Let  8  be  in  ( |e^ |,  k)  D  ( e2 ,  k),  and  define 

-A?  =  {ip-  e)  '•  ±e  €  [6,  k],  p  e  [-$,.?]} 

M  =  {(/>,e)  •'  e  €  [er,e2j,  p{e,p)  <  p  <  s } 

A0  =  {[s,s]  x  [-k,k]}\(A+  U  A~  U  A2 ) 

The  sets  .4^,  ^42,  „40  are  mutually  disjoint,  and  their  miion  is  [—A,  A]  x  [— s,  $]. 
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Theorem  5.2.  If  s  and  si  are  sufficiently  small,  then  there  is  a  6  £  ( jex ( ,  fc)  D  (e2,fc)  so 
that  the  following  hold  for  |Q|  <  s i: 

(a)  U(Q,p,  e)  >  0  if(p,  e)  €  A2.  There  exist  smooth  functions  Q\  andU~  :  AfUAo  — + 
[— Q2  and  II+  :  Af  U  Aq  — ►  [— si,si],  and  IIo  :  Ao  — ►  [— si,si],  such  that 

(b)  II (Q,p,e)  =  (Q  -  Qi(p,e))Il-(Q,p,e)  if(p,e )  £  A~ ,  where  11"  >  0;  U(Q,p,e)  = 
( Q  -  <22(p,e))n+(Q,p,e)  if  (p,  e)  £  A+,  where  II+  <  0. 

(c)  Qi(p(e),e )  =  Q2(p(e),e),e  £  [e1}  e2];  Q^p,  e)  >  Q2(p,  e),  (p,  e)  £  A0;  Q i(0,e)  = 
0,  e  £  [— k,  0];  and  Qi(p,  e)  =  0,  e  £  [0,  fc] 

(d)  II (Q,p,e)  =  (Q-  Qi(p,e))(Q  -  Q2(p,e))n0(Q,p,e)  if(p,e)  £  A0,  where  n0  >  0. 

Remark  5.3.  For  e  0  the  root  vanishing  at  p  =  0  behaves  approximately  as  —aip/a2e 
while  in  Ao  the  two  sheets  are  approximated  by  solutions  of  the  quadratic:  a4Q2  +  (ea2  + 
pa-3)Q  +pai  =  0. 

With  Theorem  5.2  in  hand,  the  analysis  of  (5.21)  is  straightforward.  Theorem  5.1 
ensures  that  either  (i)  Q'(x)  is  one-signed  on  (—00, 00),  or  (ii)  (after  a  suitable  translation) 
Q  is  even  about  x  =  0  and  monotone  on  (— 00, 0)  and  on  (0, 00, ).  Equation  (5.15)  ensures 
p  >  0;  that  is  A  <  Xd.  We  now  show  that  the  inequality  is  strict  for  a  nontrivial  solution. 
If  p  —  0  and  (p,  e)  £  Af ,  then 

O'(*)*“«3(*)n±(0(*),0,e). 

If  (i)  occurs,  then  Q(x )  approaches  a  constant  Q  as  x  — ►  00,  whence  n±(Q,0,  e)  =  0. 
However,  this  contradicts  Theorem  5.2(b).  If  (ii)  occurs,  then  Q'(0)  =  0  and  Q( 0)  ^  0, 
which  is  again  impossible.  If  p  =  0,  (p,  e)  £  Ao,  and  e  >  0,  then 

(5.25)  Q\xf  =  Q3(x)(Q(x)  -  Q1(0,e))n0(Q(x),p,e). 

For  the  case  (i),  we  have  Q(x)  — ♦  <2i(0,  e)  >  0  as  x  — *■  00,  whence  Q(x)  >0,  x  £  R,  and 
g3(x)(g(x)-g1(0,e))<0,  x€R.  Theorem  (5.2)(d)  gives  Ho  >  0,  and  this  contradicts 

(5.25) .  If  e  >  0  and  (ii)  holds,  then  g(0)  =  gi(0,e),  Q(x)  >  0,  x  €  R,  and  a  contradiction 
follows  as  before.  A  similar  argument  holds  for  e  <  0.  If  e  =  0,  then  g2(0,0)  =  0,  whence 
g'(x)2  =  g4(x)iio(g(x),o,o),  with  n0  >  0  which  is  again  impossible.  Hence,  nontrivial 
solutions  of  (5.17)-(5.20)  necessarily  have  p  >  0  or 

X  Xd, 
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so  they  are  supercritical  in  terms  of  the  speed  c  =  y/gJX. 

Next,  we  note  that  if  (p,e)  €  A2,  then  Q'(x )2  =  Q(x)2U(Q,p,e)  where  U(Q,p,e)  >  0 
by  Theorem  5.2(a).  Hence,  there  are  no  solutions  for  (p,  e)  6  A 2.  Define 

A+  =  {(0,s]x[-fc,fc]}\A2 


and 

r0  =  {(p(e)>e)  :  e  £  [e1,e2]  \  {0}}. 
which  is  a  projection  of  F  minus  the  origin. 

For  (p,  e)  £  A+,  there  is  a  ’solitary  wave’  solution  which  is  zero  at  x  =  ±oo  and  has  a 
value  at  x  =  0  equal  to  the  simple  root  of  II  nearest  zero.  When  (p,  e)  £  To,  II  has  a  double 
root  at  Qi  =  Q2  and  there  is  a  surge  solution  which  approaches  this  value  as  x  — ►  +oo. 
These  results  are  now  summarized. 

Theorem  5.4.  (a)  If  (p,  e)  £  A+,  there  is  a  unique  solution  Q(p,  e)(x)  to  (5.17)-(5.20) 
with  Q(p,e )  even  about  x  =  0.  The  functions  Q(p,e)(x)  and  Q'(p,e)(x)  have  the  same 
sign  as  e  on  (—00,00)  and  (— 00, 0),  respectively.  Furthermore  Q(p,e)(0)  equals  Q2(p,e)  if 
e  >  0  and  equals  Qi(p ,  e)  if  e  <  0. 

(b)  If  (p,  e)  £  To,  there  is  a  solution  Q(p,e)(x)  to  (5.17)-(5.20)  with  Q(p,  e)(x)  — ► 
Qi(p,e)  =  £?2(p,  e)  ^  0  as  x  — *  00.  The  functions  Q  and  Q1  have  the  same  sign  as  e  on 
(—00, 00).  Any  other  solution  of  (5.21 )  is  a  translation  of  Q(p,  e). 

The  solutions  in  (a)  give  rise  to  solitary  wave  solutions  u;(x,  y )  while  those  in  (b)  give 
surges.  The  solutions  in  Theorem  5.4  can  be  analyzed  with  a  scaling  argument.  For  the 
case  (p,  e)  €  T0, 

Q\xj>  =  Q\x)(Q(x)  -Q10>,(0)Jn0(Q(*),j>,<0. 

Let  Q\(p(e),e)  =  q  and  define  F(x)  by  the  relation  Q(p,e)(x)  —  qY(\q\x)  whence 

Y\xf  =  r2(i)(K(x)  -  i)2n„(?r,p,e). 

Note  that  Y'{x)  >  0  and  F(x)  ->  1  as  i  ->  00.  As  e  — »  0,  <7  — ►  0,  and  so  Y(x)  will  behave 
like  a  solution  to 

y'{xf  =  y2(x)(y(x)  -  l)2IIo(0,0,0) 

=  y20)(y0)  -  i)2a4(o,o). 
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Using  (5.18)  one  has 


where 


!/(*)  =  i(l+tanh(^)), 


a\p>  fll/4h  x«3/2V 


pl/4(l  +  p3/2)' 

For  small  e,  Qi(p,  e)  is  given  by  (5.23),  whence 


/-  w  \  ^  V3  y/p\/l  +  p3^  f,  ,  x  ^r>/3  P1/4 _ jejx  . 

(5.  6)  Q(p,e)(x)  12  (1  +  v^9/2  4  (1  +  y^)3/2  ^Tf?/'2  ' 

Note  by  (5.23)  that  |e|  =  0(^/p)  and,  using  (5.16)),  R(x,y)  =  0(|p||Q|  +  \Q\ 2)  =  0(e2). 
Hence,  (5.26)  gives  an  expansion  for  the  bore  solution  w(x ,  y)  =  Q(x)t(y)+R(x ,  y)  accurate 
to  order  e,  or  equivalently,  order  \/\d  —  A. 

Suppose  (p,  e)  €  A+,  and  consider  the  behavior  of  Q(p,e)  as  p  — ►  0  for  /ixed  e.  Now 
one  sets  Q(p,  e)(x)  =  tfY(-\/igjx),  where  q  =  Qi(p,  e)  for  e  <  0  and  y  =  Q2(p,e)  for  e  >  0, 
and  deduces  that  Y  behaves  like  a  solution  to 


y'(x)2  =  m3h3\e\y2(x)(l  -  y(x)), 


that  is, 


It  follows  that 


,(x)  = 

«(P,e)(x)^^sech^(<^Mk). 


emh  2 

Note  that  this  solution  tends  to  zero  like  |p|  =  jAj  —  A|  for  fixed  e,  while  the  bore  solution 
in  (5.27)  has  the  slower  decay  0{y/\p\). 

For  a  fixed  (p,  e)  €  A+,  Theorem  5.4  gives  the  existence  of  a  homoclinic  orbit  Q(p,  e ) 
to  (5.17)-(5.20)  with  Q(p ,  e)(0)  =  Qi(p,  e).  This  is  a  closed  loop  in  the  ( Q ,  Q')  phase-plane. 
One  can  show  easily  that  the  trajectory  encircles  only  one  critical  point  ( Q,Q )  (that  is,  a 
non-zero  constant  for  which  the  right-hand  side  of  (5.17)  vanishes).  It  is  then  immediate 
that  the  interior  of  the  trajectory  is  filled  with  small  periodic  solutions  of  (5.17)  which  give 
rise  to  periodic  (in  x)  solutions  w(x,y). 
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